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We study the DVCS amplitude within the color dipole approach. The light-cone wave function 
of a real photon is evaluated in the instanton vacuum model. Our parameter free calculations are 
able to describe HI data, both the absolute values and the t-dependences, at medium-high values of 
. The dependence is found to be sensitive to the choice of the phenomenological cross section 
fitted to DIS data. 
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■ I. INTRODUCTION 

CN ■ 

^ ' During the last decade hard exclusive reactions, such as deeply virtual cornpton scattering; (DVCS), 7* + P ~^ 7 + V, 
, have been a subject of intensive theoretical and experimental investigation [E 0, S 0, H, H, 0, B B EO) El) [3 El, [iJ, [3l ■ 
Particular interest deserves the Bjorken kinematics, 

: - 9' = W^Q^ » KlcD, W' = {p + qf » AlcD, (1) 

XB = const, t = A^ = {p'-pf^Q\ (2) 

I ' where q is the momentum of the virtual photon, p is the momentum of the inital hadron, p' is the momentum of the 
final hadron, and t is the momentum transfer. In this kinematics the DVCS amplitude is factorized V, "s] into the 
(-H ' convolution of a perturbative coefficient function with a soft matrix element-generalized parton distribution (GPD) 
I— I, of the parton inside the target. 

While for the region a; 1 the dominant contribution to the DVCS amplitude comes from the generalized parton 
distributions (GPD) for quarks, we know from experience with deep-inelastic scattering (DIS) that in the small- a: 
kinematics (large cm. energy W) the density of partons (especially gluons) increases and the dominant contribution 
comes from the gluon sector. 

^\ . In this paper we employ the color dipole approch introduced in [l^. \v^. The central object of the model is the 
dipole scattering amplitude, A{x,r). For DIS, only the imaginary part of the amplitude contributes, and the result 
, becomes especially simple [1^, [l^ , 

g ■ a^.p{x,Q'') ^ J dzj d\ \^{z,r,Q^)\''agg{x,r), (3) 
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-the DIS cross-section equals the dipole cross-section aqq{x,r) "averaged" with weight j^" (z, r, Q^) | , where 
^ (z,r, Q^) is the virtual photon light-cone wave function, which depends on transverse separation r and fraction 
5^ z of the light-cone momentum carried by the quark (see Sectio n IIIII for rigorous definition). 

Generahzation of ^ to the DVCS case is straightforward although one has to pay attention to some 

subtle points. First, one has to deal with the nonperturbative contribution of large-size dipoles. While it also exists 
in the forward case (DIS), the problem is exacerbated in the off-foward case (DVCS) due to presence of the real 
photon in the final state. A priori, one cannot say how important this contribution is in a p articular process. As 
a first approximation, one can just ignore the nonperturbative effects, as was done in [l^, [20, [U, [l^, and evaluate 
everything perturbatively. While for DIS this approximation is partially justified due to high virtualities in both 
vertices, for DVCS its validity is questionable, since the final photon is real. In the present evaluations we rely on the 
instanton vacuum model for the photon wave function. 

Another complication is that in contrast to DIS, the DVCS amplitude is sensitive to both the real and imaginary 
parts of the dipole amplitude. In this evaluations we restore the real part using the prescription proposed in [23| (see 
Section [n] for more details). 
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The paper is organized as follows. In Section |TT] we review the color dipole model and present the parameters used 
in numerical evaluations. In Section Hill we evaluate the photon wave function within the instanton vacuum model. In 
Sections HV] and |V] we present results for the DVCS cross-sections, compare with available data, and draw conclusions. 



II. COLOR DIPOLE MODEL 



As was mentioned above, the color dipole model is valid only in the region of small-x, where the dominant contri- 
bution to the DVCS amplitude comes from gluonic exchanges. The general expression for the DVCS amplitude in the 
color dipole model has the form. 



A- {xB,t, Q2) oc e«e(/) J d(3id(32d\id?r2^^'^ (/^z, ^2, = O) A% (/3i, rl; /32, fa; Q^ A) ^^^^ (/3i, Q') , (4) 

where /3i^2 are the light-cone fractional momenta of the quark, and fi.2 are the transverse distances in the final and 
initial states respectively (in off-forward kinematics they are no longer equal), A is the momentum transfer in DVCS, 
■4'^(...) is the scattering amplitude for the dipole state, and indices refer to polarizations of inital and final 

photons. In general case the amplitude A'^{...) is a nonperturbative object, with asymptotic behaviour for small r 
determined from the pQCD [iGj: 

A\...)^r^ (5) 

(up to slowly varying corrections ln(r)). 

In our approach we use a model [24|, [2^ HI] for the partial amplitude of the scattering of the color dipole on the 
nucleon, which has the form 



^f, (/3i,ri;/32,r2;Q2,A) « 5 {[3^ - P^) 5 {f^ - r2) d^be'^hm fg{r,X Pi) 



(6) 



Imjg{r,h,f3) - — 



d^kd^A 



^ ^ ~j ( ^ "2" 



- s 2 



[x,k, a) 



ib-A 



'iPr-(k- 



i{l-f3)r-[k- 



iPr-(k+i 



-i(l-/3)f- fe+4 



(7) 



where b is impact parameter; ^ T'^^" t) ^^'■^ transverse momenta of the incoming and outgoing quarks 
(x, k, Aj is the unintegrated GPD for gluons [i^ . 

(x, k, 



A;2 



H \^x,k,A 



(8) 



dz 



2' 2 



7+ -Co 



2' 2'2'2 



2' 2 



where £00 (x; ij] is the Wilson factor required by gauge covariance. For this GPD we use the gaussian parameteriza- 
tion [li,!!!,!!^. 



(x, k, a) 



Rlix) 



A2 



exp 



-B{x)A' 



(9) 



where ao{x), Rq{x), B{x) are free parameters fixed from the DIS and scattering data. We shall discuss them in 
more detail in Section IIVI The parameterization ([9]) does not depend on the longitudinal momentum transfer and 
decreases exponentially as a function of A^. Since it is an effective parameterization valid only in the small-x region, 
we do not assume that it satisfies general requirements, such as positivity 27 1 and polynomiality constraints. 

The prefactor guarantees convergence of the integrals in equation (j6]) . In the forward 

limit, the amplitude ^ reduces to the saturated parameterization of the dipole amplitude proposed by Golec-Biernat 
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FIG. 1: Comparison of fc-integrated forward limit qbk {x,Q^) of parameterization with realistic parameterizations: 
GRV HI], MRST [H and CTEQ5 [H. 

and WiisthofF (GBW) [1^. The correspondin g fc- integrated forward limit of the GPD ^ is compared with realistic 
parameterizations of the gluon PDFs - GRV [H, MRST [30| and CTEQ5 [Hi in the Figm-e We find that in 
the region x < 10~^ the parameterization ([9]) agrees, within experimental and theoretical uncertainties, with these 
realistic parameterizations, but deviates for larger x. 
Taking the integrals in ([5]), we arrive at 



Im Jg{r, A, /3) = / <eh e^'^Im !g{r, b, P) = 



fio(^) 



Zol£l exo f ^ , 

qq \' T^Tl^J 4 '=^P \^ \^ 2 ' 16 



(10) 



Evaluation of the real part of the amplitude is quite straightforward. As it was discussed in [23|, if the limit 
lims^oo ( "^™"^ ) exists and is finite, then the real part of the amplitude is related with imaginary part by 



7^e ^ = s" tan 



n f ^ d 

2 r"'+9h^ 



ImA 



(11) 



As it will be shown in Section IIVI the imaginary part in the color dipole model indeed has a power dependence on 
energy, ImA{s) ~ s". The formula pTjl in this case may be simplified to 



UeA 
ImA 

Then for the DVCS amplitude we finally obtain 



— tan (^(a — 1)) := -q. 



(12) 



A^ « (77 + z)e«(g')eL^'(9) / d\ / r, = 0)vI/W(/3, r, Q^) /^(^- A, /5) 



(13) 

where in psp the upper indices (ij) refer to polarizations of the final and initial states, and we evaluated the real 
part of the DVCS amplitude according to (llip . For the cross-section of unpolarised DVCS amplitude, from (fT3|) we 
obtainii 



da 

Itt 



167r 



1+^ 
167r 



E 



(14) 



Since the helicity-flip scattering is a higher-twist effect, it is suppressed in the large-Q^ kinemtics. Thus we may 
further simplify (|14p to 



4 



da 
lit 



1 + ?7 



16 



-E 



0)*«(/3,r,Q2),/^/^^(^,A,/3) 



(15) 



Formula |T5|) is the final result which will be used for numerical evaluation of the DVCS amplitude. It is important 
to note that in the small-r region the amplitude (fTO|) behaves as ^ ~ r^, therefore the corresponding cross-section (fT4|) 
is falling as 1/Q^ at large-Q^, in agreement with the general analysis P, 



III. PHOTON WAVE FUNCTION IN THE INSTANTON VACUUM 



In this section we would like to provide some details of evaluation of the wave function in the instanton vacuum model 
(IVM) (see [13, [s^, [1^ [IHl and references therein). In the leading order in Nc^ the model has the same Feynman 
rules as in the perturbative theory, but with momentum-dependent quark mass ^{p) in the quark propagator [33l | 



S{p) 



1 



p - /i(p) + iO' 



and nonlocal interaction vertex of vector current |35l | 

V = -^V = v-M (G^(p, q)f{p + q) + G,,{p + q, -q)f{p)) v^^iq), 



(16) 

(17) 
(18) 



n=0 



where p, p + q are the momenta of the incoming and outgoing quarks respectively, and f,i_ii...i_i„{p) = 
f{p) / dp^-^...dp^^. Using symmetry properties of the the last term in (|17p and properties of the fimction G^(p, q), 
it is possible to cast (|17l) to the equivalent form 



V ^v^'{q)\^^,~{2p^ + q^) 



which is frequently used [371, |38j . 

The mass of the constituent quark has a form [SS 



M{P{p + q)-p{p)y 

ip + qY - p^ 



(19) 



M(p) = m-fM/2(p), 



(20) 



where m « 5 MeV is the current quark mass, and M « 350 MeV is the contribution of the instanton-induced effects. 
In the limit p — > cx) the instanton-induced nonlinear formfactor f{p) falls off as ~ p-, so for large p 3> where 
p « (600A/ey)~^ is the average instanton size, the mass of the quark ^(p) k, ra and the vector current interaction 
vertex V k, i). However, we would like to emphasize that all the correlators get contributions from both the soft and 
the hard parts, so even in the large-Q limit the instanton vacuum wave function is different from the well-known QED 
result. 

The overlap of the initial and final photon wave functions in (llSp was evaluated according to 



«'«*(/3, r, Q2 = iP, r, Q') = ^ (/3, r*, 0) Ir (/3, r, Q') , 



(21) 



where the summation is done over possible polcLriza.tion sta-tes T^em — {o'/^; T/^TS; ^/^^^ 

}, and Ir corresponds to one of 

the matrix elements 



dz 
Yn 

dz~ 



2n 



2' 2r''^^^V2'2 



7(9) 
7(9) 

7('Z)\ • 



(22) 
(23) 
(24) 
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FIG. 2: In the leading order there is only one diagram which contributes to the wave function (|25[). 



Notice also that in the final state in [21] we should use r* — + 



g±-r± 
1+ 



9+ 



which takes into account 



that the final photon has 7^ whereas the components of the wave function (|22II24[) are defined in the reference 
frame with q± = 0. 

The wave functions corresponding to matrix elements (I22II24|) were evaluated in [s^. In the leading order in Nc 
one can easily obtain for the components /r 



Ir = 



(^S{p)VS{p + q)Te 



(25) 



which corresponds to the diagram shown in the Figure ^ . The evaluation of (PS)) is quite straightforward and in the 
reference frame with q± = 0, ^x^l) ~ ^\ (?) yields 



(26) 



+ g)" + P'jP + <l)^^ - {P^ +P-<1- Kp)t4P + g)) 

^ {2l3q+p- -pl-^I^ (p) + iO) (2 (1 - P) q+ {p- -pl-^'^{p + q) + jo) 

_ M {G, {p, q)f{p + q) + f jp) {p + q, -q)) {fi{p){p + q)^ + fi{p + q)pf,) \ 
{2Pq+p- -pl^fi^ (p) + iO) (2 (1 - /?) q+ (p- - 1^) - pI -fiHp + q) + lO) ) ' 



d^p 
(2^ 



-f})q+{ 


P - 






-ir± (p+f 


)x 


Pj 





z r 
2' 2 



7'(g) 



(27) 



(2/59+p_ - pI ~ + iO) (2 (1 - /3) q+ {p- - 1^) - Pi - m'(p + g) + *0) 



7'(g) 



(28) 



--.4^'(g)/|f.e-(-^)x 

(gi^g^A - g^gt^A) + (a^(p + g) - m(p)) (p^g,.A " Pt/ff^A) + M (Ga (p, q)f{p + q) + f jp) Gx {p + g, -g)) (p^g,. - Pt/g^) 

(2/3 g+p- - pI - (p) + lO) (2 (1 _ ;3) (7+ (p- ~ 1^) - pi ~ f^Hp + q) + io) 
We evaluated (|26ll28p numerically. 



IV. RESULTS 



Here we present the results of the evaluation of the DVCS differential cross section with different models for the 
partial dipole amplitude (r , b) , and we also use a photon wave function calculated either perturbatively, or within 
the instanton vacuum model. First we test the amplitude based on the GBW model. We also perform calculations 
with the parameterization proposed by Kowalski and Teaney (KT) for the saturated dipole cross section 13911 . In 
addition, we try an energy dependent KST parameterization for the dipole cross section, proposed in [25l. l26l. lioj. 
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A. GBW based partial dipole amplitude 



First of all, we made evaluations in the GBW model extended to the 6-dependent partial amplitude f^{r,b). 
This parameterization is fitted to DIS data at large and small Bjorken x. The parameters in Eq. ([9]) read, 
ao{x) = 23.03 mb^const, Ra{x) = 0.4 x (a;/a;o)'' "'*fm. where xq = 3.04 x 10"''. The parameter B{x) in Eq. (HI), is 
related to the t-slope of the differential cross section of highly virtual photoproduction of vector mesons 24, 25, 26], 

Bix)=B^,p^pp^^Rl{x) (29) 

We use the e xpe rimental value of the slope Bj-p^pp 

Following [28|, we use the perturbative wave function of the photon with constituent quark mass 140 MeV. From 
the left panel of Figure [3] we see that in the small- a;^ region the cross-section is proportional to a power of xg , 
da/dt ~ Xg, where the power a was obtained by fitting the a;^— dependence in the range xb G (10^^,10"'^) and 
slowly depends on (Q^, t). 

In the right panel of Figure [3] we compare the model with data from the HI experiment at HERA . Although 
the model gives a reasonable description of the data at moderate values of Q^, the discrepancy increases at higher Q^. 
As it was discussed in SectionlTTl the modeled DVCS cross-section falls as l/Q^, while the data behave approximately 
as 1/Q3. 




FIG. 3: Xb— and t— dependence of the DVCS cross-section in GBW parameterization 



In order to identify the source of the disagreement, we replace the perturbative photon wave function by the one 
calculated in Section IIIII within the model of instanton vacuum. The results are depicted in Figure 21 We see that 
the cross sections do not change very much compared to the previous calculation. This model gives a reasonable 
description of the cross-section for moderate Q^, however grossly underestimates the data at high Q^. 

While the GBW dipole cross section does a pretty good job describing DIS data, and electroproduction of vector 
mesons 43, 44], it fails to explain the observed dependence in DVCS. This fact shows that DVCS provides a rather 
sensitive test for models. 
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da/dt , nb/GeV^. t=-0.1 GeV^ 



da/dt , nh/GeV^, W=82 GeV 



Q^=4 GeV^ da/dt ~ Xg^^f 





10-' - 



10-^ 



10-* =- 



■ Q-=8 GeV^ 
. Q-=15.5 GeV^ 
Q^=25 GeV' 



I I I I I I I I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Xb -t, GeV' 

FIG. 4: xb— and t-dependence of the DVCS cross-section in GBW parameterization, with realistic photon wave function. 

B. KT parameterization 

Another form of impact parameter dependent partial dipole amplitude, which has correct behavior at small r and 
the saturated shape, was proposed in [39|, 



Im/(r, 6) = 2 1 - exp 



(30) 



where the scale n = 0.77 GeV + 4/r , and the gluon distribution function g ) was fitted to DIS data. The 
nucleon profile function in [soj, Tj^^b), has a simple form 



TN{b) 



1 



2nBn 



(31) 



where the slope parameters Bg — AGeV~^ was fitted to data on electroproduction of vector mesons. 

In the small-r^ limit this function corresponds to ordinary gluon PDF g (a;, /i^) . However, the latter may be different 
from the results of DGLAP analyses of data, since Eq. (|30p is supposed to include the saturation effects. On the 
other hand, the eikonalization used in (j30|) is quite a rough procedure at large r where saturation is at work. A more 
accurate formula should include a convolution of the dipole amplitude with T/v(&), rather than the simple product. 
This is why Eq. ([50)1 misses correlations between b and f , which are present in ([7]) . 

Results for the DVCS differential cross section, calculated with this parameterization and the realistic nonpertur- 
bative photon wave function, are presented in Figure [51 
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C. Energy dependent (KST) parameterization 

While in the large-Q^ kinematics scahng implies that the dipole amplitude is a function of Bjorken x, for smaller 
(soft photons) scaling does not work, and for real photons Bjorken x is not an appropriate variable. For this 
kinematics the photon energy s should be used instead of xb- An s-dependent parameterization with saturated form 
analogous to the GBW was proposed in [40| , for the description of the real photo production and absorption, and for 
DIS at small . One should replace all the a;-dependent functions in Eq. ([9]) by s-dependent ones, with 

a,{s)^aZ{s)[l + \^y (32) 

where a'l^^{s) = (Sq + Si ln^(s/so)) with Sq 20.9 to6. Si = 0.31 m6, and sq = 28.9 GeF^^ ^^le total pion- 
proton cross section, and {r'^h) ~ QA4:fm? is the pion charge radius. Correspondingly, Ro{s) — 0.88(5/51)"-^^ /m; 
si = 1000 GeF^. In this case the parameter B{s) is related to the t-slope of elastic up scattering [H, HH], 

B{s)^B:i'{s)-\Rl{s)-\{rl,). (33) 

Naturally, we use the nonperturbative photon wave function Eq. (j2ip . 

The numerical results for the DVCS cross section are depicted in Figure [51 Comparison with data plotted in the 
right panel show that this model also leads to a too steep dependence of the cross section, although its absolute 
value at = 8 GeV^ agrees with the data. 
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dc/dt , nb/GeV^. t=-0.1 GeV^ da/dt , nh/GeV^, W=82 GeV 




10" 10^ 10^ 10' 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Xb -t, GeV^ 

FIG. 6: xb— and f-dependence of the DVCS cross-section, in KST parameterization, with realistic photon wave function. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we evaluated the DVCS cross-section relying on a nonperturbative wave function of the real photon 
calculated in the instanton vacuum model. The ratio of DVCS cross sections calculated with the IVM and perturbative 
wave functions is shown in the left panel of Figure [7] as function of t for different values of Q^. 




-t, GeV^ Q\ GeV^ 



FIG. 7: Left: Ratio of the DVCS cross-section with pQCD and IVM wave functions in GBW model. One can clearly see that 
for Q2 ~ 8GeV^ the difference might reach up to 50fimodels (GBW and KST). 

This comparison demonstrates a considerable modification of the cross section up to factor two, dependent on 
kinematics. For small-Q^ ~ 1 GeV^ one should use instanton wave functions both for the initial and final photons. 
Note that even for ~ 1 GeV^ both cross-sections dajvAi/dt and dapQco/dt decrease quite fast as a function of t, 
approximately as e^*, as one can see from the previous Figures [3][5] The difference is due to the higher- twist effects 
which are amplified in the small-Q^ region and become pronounced in the ratio of the cross-sections. 
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We also tested several models for the dipole cross section. Although all the parameterizations under discussion 
have been fitted to DIS data from HERA, not all of them are successful in describing the DVCS cross section, 
especially the observed dependence. In particular the dipole partial amplitude based on the popular GBW 
parameterization of the total cross section leads to a too steep dependence of the DVCS cross section, calculated 
with both IVM and perturbative photon wave functions. In this model the cross-section decreases in the measured 
interval Of = 8 — 25 GeV^ as l/Q^, like one should expect at very large in accordance with the general large-Q^ 
analysis [l|, i, 0], whereas HI data exhibit an approximate ~ l/Q'^-behaviour. We compare the Q^-dependences 
calculated with different models for the photon wave function and dipole partial amplitude, at fixed xb and t in the 
right panel of Figure [T] We also show in the Figure the results of the fit to the calculated cross sections within the 
interval measured in the HI experiment. We conclude that although the absolute value of the cross-section is 
sensitive to the nonperturbative effects (large-size dipoles), the Q^-dependence does not vary much. 

Apparently, one can achieve a weaker dependence by replacing the standard quadratic behavior of the 
dipole cross section predicted by pQCD at small r [l^ by a smaller power of r. Such a model [11] was considered 
in to describe the Q^-dependence of HI data for DVCS cross section by modifying the small-r behaviour to 
Ad ^ rpj^jg cxplaius wliy the model gives reasonable description for the Q^-dependence of the cross- 

section. 

Remarkably, the KT model [l^l, considered above in section |IVB] provides a reasonable description of the measured 
Q^, although it has the same small-r behaviour Ad ~ r'^ as GBW. This happens because the higher- twist effects in 
this model arc more pronounced than in GBW. Straightforward expansion of the dipole amplitude Ad in the KT 
parameterization yields (we consider the case Aj_ = for the sake of simplicity). 

Ad -r^ [l~XKTir)r^+0{r^]), (34) 

where the constant XKri^) depends on r only logarithmically, and for typical (r) ~ (Q)^^ ~ 0.1 GeV^"^ we have 
Xxrir) ^ 0.159 GeV^, which is larger than in GBW. 

We would like to emphasize that although it is possible to fit the DIS data with a parameterization containing 
sufficient number of free parameters, not all of the available parameterizations are able to describe the DVCS data, 
and a systematic study of the higher-twist corrections is necessary. In this paper we addressed one of the possible 
sources-the higher-twist corrections to the wave function of the initial photon and considered realistic wave function 
for the final (real) photon. Other sources of higher-twist corrections include contributions of higher-twist quark-gluon 
operators, such as (p' \qq\p} , (p' |qG"G+c((7| p) etc. Contrary to DIS and heavy meson electroproduction, in most of 
other processes the large-r behaviour of the amplitude is important and gives contribution comparable to the small-r 
region. That's why in modelling of the large-r behaviour of the amplitude we should refer to some microscopic model 
rather than to a scmiphcnomcnological approach. 
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